Proton-neutron interaction near closed shells 
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Odd-odd nuclei around double shell closures are a direct source of information on the proton- 
neutron interaction between valence nucleons. We have performed shell-model calculations for dou- 
bly odd nuclei close to 208 Pb, 132 Sn and 100 Sn using realistic effective interactions derived from 
the CD-Bonn nucleon-nucleon potential. The calculated results are compared with the available 
experimental data, attention being focused on particle-hole and particle-particle multiplets. While 
a good agreement is obtained for all the nuclei considered, a detailed analysis of the matrix elements 
of the effective interaction shows that a stronger core-polarization contribution seems to be needed 
in the particle-particle case. 

PACS numbers: 21.60.Gs, 21.30.Fe, 27.60.+j, 27.80.+W 



I. INTRODUCTION 

Odd-odd nuclei around double shell closures provide 
the best testing ground for the matrix elements of the 
proton-neutron interaction between valence nucleons. In 
this context, of special interest are nuclei in the close 
vicinity to 208 Pb, 132 Sn and 100 Sn, which show strong 
shell closures for both protons and neutrons. From the 
experimental point of view, the 208 Pb neighbors have 
been extensively investigated and a large amount of data 
is available for them. On the other hand, the 132 Sn and 
100 Sn neighbors lie well away from the valley of stability, 
which makes it very difficult to obtain information on 
their spectroscopic properties. In recent years, however, 
substantial progress in the development of spectroscopic 
techniques has opened the way to the exploration of nu- 
clear structure in the regions of shell closures off stability. 
This has led to new experimental data for the immedi- 
ate odd-odd neighbors of both 132 Sn and 100 Sn. The 
information presently available for nuclei of this kind, 
although still rather scanty, is of great value for the un- 
derstanding of the effective proton-neutron interaction 
in these two regions as well as for a comparison with the 
proton-neutron interaction in the 208 Pb region. 

During the past several years, we have studied a num- 
ber of nuclei around 208 Pb, 132 Sn and 100 Sn within the 
framework of the shell model employing realistic effective 
interactions derived from modern nucleon-nucleon (NN) 
potentials. In most of these studies, however, we have 
been concerned with nuclei having few identical valence 
particles or holes. A summary of the results obtained 
through 2000 is given in Ref. [n. 

More recently 0, , we have turned our attention to 
nuclei with unlike valence nucleons to try to gain infor- 
mation on the neutron-proton effective interaction. The 
main aim of this paper is to report on some selected 
results of our current work along these lines, which we 
have obtained starting from the CD-Bonn free NN po- 
tential 4]. In particular, we consider the six odd-odd 
nuclei 208 Bi, 210 Bi, 132 Sb, 134 Sb, 98 Ag and 102 In, and 
focus attention on the particle-hole and particle-particle 



multiplets, which offer the opportunity to test directly 
the matrix elements of our calculated realistic effective 
interaction. 

To place this study in its proper perspective, it should 
be mentioned that the particle-hole multiplets in the 
208 Pb region were the subject of great theoretical and 
experimental interest 0, 0, @ some thirty-five years 
ago. In particular, the spectrum of 208 Bi was extensively 
studied through pick-up and stripping reactions and sev- 
eral proton- neutron hole multiplets were identified 0, |(| . 
The pattern exhibited by the observed multiplets was 
well reproduced by the shell-model calculations of Ref. 
, where particle- hole matrix elements derived from the 
Hamada- Johnston NN potential were used. 

Despite these early achievements in the understanding 
of the effective proton-neutron interaction around closed 
shells, little work along the same lines has been done 
during the last three decades. The new data which are 
becoming available in the 132 Sn and 100 Sn regions, as well 
as the prospect of spectroscopic studies of unstable nuclei 
opened up by the development of radioactive ion beams, 
has motivated us to undertake the present study and per- 
form shell-model calculations making use of a modern 
NN potential and improved many-body methods for de- 
riving the effective interaction. 

The outline of the paper is as follows. In Sec. II we 
give a brief description of our calculations. Our results 
are presented and compared with the experimental data 
in Sec. Ill, where we also discuss the effects of the core 
polarization. In Sec. IV we draw the conclusions of our 
study. 



II. CALCULATIONS 



In our calculations for the Bi isotopes we assume that 
208 Pb is a closed core and let the valence proton and neu- 
tron hole occupy the six single-particle levels 0/ig/2, l/7/2> 
0«i3/2, 1/5/2, 2p 3 / 2 and 2pi/ 2 of the 82-126 shell, while for 
the valence neutron the model space includes the seven 
orbits I39/2, Oin/2, 0ji5/2, 2d 5 / 2 , 3si/ 2 , I57/2, and 2d 3/2 - 
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Similarly, for the Sb isotopes we assume that the valence 
proton and neutron hole occupy the five single-particle 
levels O57/2, lds/2, 1^3/2, 2si/a and Ohu/2 of the 50-82 
shell, while the valence neutron outside the 132 Sn core 
occupies the orbits of the 82-126 shell. As regards 98 Ag 
and 102 In, we assume that the valence neutrons outside 
the 100 Sn core occupy the orbits of the 50-82 shell, while 
for the proton holes the model space includes the four 
orbits O59/2, lpi/2, IP3/2 an( i O/5/2 of the 28-50 shell. 

As regards the choice of the single-particle and single- 
hole energies, we have proceeded as follows. For the Bi 
isotopes we have taken them from the experimental spec- 
tra [l3 of 209 Bi, 209 Pb and 207 Pb. In the same way, for 
the Sb isotopes we have made use of the experimental 
spectra of 133 Sb [HI, 133 Sn [j3 and 131 Sn In the 

spectra of the two former nuclei, however, some single- 
particle levels are still missing. More precisely, this is 
the case of the proton 2si/2 and neutron 0?i3 /9 levels, 
whose energies have been taken from Refs. and [l4| . 
respectively, to which we refer for details. As regards the 
neutron hole energy e^ 11/2 in 131 Sn, we have adopted the 
value of 0.100 MeV which has been recently suggested in 
Ref. [13. This is somewhat smaller than that reported 
in Ref. - For 98 Ag and 102 In the single-particle and 
single- hole energies cannot be taken from experiment, 
since no spectroscopic data are yet available for 101 Sn 
and 99 In. Therefore, we have taken them from Ref. 0] 
and Ref. where they have been determined by an 

analysis of the low-energy spectra of the Sn isotopes with 
A < 111 and of the N=50 isotones with A > 89. For com- 
pleteness, our adopted values for the single-particle and 
single-hole energies are reported in Tables I-III. 

As already mentioned in the Introduction, in our shell- 
model calculations we have made use of a realistic effec- 
tive interaction derived from the CD-Bonn free nucleon- 
nucleon potential Q|. This high-quality NN potential, 
which is based upon meson exchange, fits very accurately 
(x 2 / datum ~ 1) tne world NN data below 350 MeV 
available in the year 2000. 

The shell-model effective interaction V e g is defined, as 
usual, in the following way. In principle, one should solve 
a nuclear many-body Schrodinger equation of the form 



HV, = EM 



(1) 



with H = T + Vat at, where T denotes the kinetic en- 
ergy. This full-space many-body problem is reduced to a 
smaller model-space problem of the form 



PH eS P^i = P(H + V eS )P^i = E.P^f, 



(2) 



Here Hq = T + U is the unperturbed Hamiltonian, U be- 
ing an auxiliary potential introduced to define a conve- 
nient single-particle basis, and P denotes the projection 
operator onto the chosen model space, 



d being the dimension of the model space and the 
eigenfunctions of Hq. The effective interaction V B s oper- 
ates only within the model space P. In operator form it 
can be schematically written as 



(3) 
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V eS = Q-Q' J Q+Q' J Q J Q-Q' J Q J Q J Q+ ... , 

where Q, usually referred to as the Q-box, is a vertex 
function composed of irreducible linked diagrams, and 
the integral sign represents a generalized folding opera- 
tion. Q' is obtained from Q by removing terms of first 
order in the interaction. Once the Q-box is calculated, 
the folded-diagram series of eq. (4) can be summed up 
to all orders by iteration methods, as, for instance, the 
Lee-Suzuki one [Hill. 

A main difficulty encountered in the derivation of V c s 
from any modern NN potential is the existence of a 
strong repulsive core which prevents its direct use in 
nuclear structure calculations. This difficulty is usually 
overcome by resorting to the well known Brueckner G- 
matrix method. Here, we have made use of a new ap- 
proach I2IJ w hich has proved to be an advantageous al- 
ternative^], I22I |23| to the use of the above method. 
The basic idea underlying this approach is to construct 
a low- momentum NN potential, Vi w-k, that preserves 
the physics of the original potential Vnn up to a cer- 
tain cut-off momentum A. In particular, the scattering 
phase shifts and deuteron binding energy calculated by 
Vnn a re reproduced by Mow-k- The latter is a smooth 
potential that can be used directly as input for the cal- 
culation of shell-model effective interactions. A detailed 
description of our derivation of Vi ow -k can be found in 
Ref. |2l|, where a criterion for the choice of the cut-off 
parameter A is also given. We have used here the value 
A = 2.1 fin- 1 . 

Once the Viow-k is obtained, the calculation of the ma- 
trix elements of the effective interaction is carried out 
within the framework of the folded-diagram method out- 
lined above. The key feature of this approach is that 
there is no need for G-matrix procedures to eliminate the 
effects of a strongly repulsive core, since the Vj ow -k can 
be used directly in the calculation of the vertex function 
Q-box. The calculation of Q, which is in principle an 
infinite sum of irreducible diagrams, can only be made 
approximately by selecting certain classes of diagrams. 
In our calculations we have included in Q all the 1-body 
and 2-body diagrams up to second order in Viow-k- For 
instance, for the particle-particle case these are precisely 
the diagrams shown in Ref. [24|. In this connection, it 
should be pointed out that the proton-neutron matrix el- 
ements for 208 Bi, 132 Sb, 98 Ag and 102 In have been explic- 
itly derived in the particle-hole formalism. A description 
of the derivation of the particle- hole effective interaction 
is given in Ref. |2( ■ 

To summarize, there are three main steps in our deriva- 
tion of Veff. We first derive the low- momentum Vi ow _k 
from the CD-Bonn potential, then calculate the Q-box 
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including diagrams up to second order in V[ ow _k, and fi- 
nally obtain V e g by summing up the folded diagram series 
(4) by means of the Lee-Suzuki iteration method. 



III. RESULTS AND COMPARISON WITH 
EXPERIMENT 

In this section, the results of our calculations for 
the proton-neutron multiplets around 208 Pb, 132 Sn, and 
100 Sn are presented and compared with the experimen- 
tal data. In the first two mass regions, both proton- 
neutron hole and proton-neutron multiplets are consid- 
ered by studying 2 os,2io Bi and I32,i34g b These nuclei, all 
having only two valence nucleons, are the most appropri- 
ate ones to study the effects of the proton-neutron inter- 
action. As regards doubly magic 100 Sn, the two immedi- 
ate odd-odd neighbors for which some experimental in- 
formation is presently available are 98 Ag and 102 In, both 
with four valence nucleons. We have therefore considered 
these two nuclei to study the proton hole-neutron mul- 
tiplets in this region. We present our results separately 
for each mass region in the three following subsections. 
All calculations have been performed using the OXBASH 
shell-model code HH. 



A. 208 Pb region: 208 Bi and 210 Bi 

The particle-hole multiplets in 208 Bi were long ago 
found to exhibit a peculiar behavior 0, |(| , namely the 
states with the minimum and maximum J have the high- 
est excitation energy, while the state with next to high- 
est J is the lowest one. As mentioned in the Introduc- 
tion, the results of the shell-model calculations of Ref. 
turned out to be consistent with this pattern. 

In Fig. I our results for three multiplets are reported 
and compared with the experimental data [ill ] . Note that 
the energies are relative to the 5 + ground state, which is 
a member of the doublet irhg/^^Pi/^- We see that the 
calculated energies account for the pattern of the exper- 
imental multiplets, the quantitative agreement being re- 
markably good with discrepancies well below 100 keV for 
almost all the states. It should be pointed out that our 
calculations give a better description of the experimental 
data than the early realistic calculations of Ref. |t| • 

Let us now come to 210 Bi. In Fig. 2 the cal- 
culated TThg/ 2 i / gg/2 multiplet (represented by open cir- 
cles) is shown and compared with the experimental data 
[Toj . This multiplet shows a breakdown of the Nordheim 
strong rule [2(j in that the l - state is the ground state 
while the 0~ one lies at about 50 keV excitation energy. 
We see that our calculation predicts for the ground state 
J OT = 0~ with the 1~ state at about 150 keV excita- 
tion energy. This is, however, the most significant dis- 
crepancy, all other experimental excitation energies being 
reproduced within few tens of keV. As a whole, the pat- 
tern of this multiplet as compared to that of the particle- 



hole multiplets is still parabolic, but concave downward. 
Of course this stems from the fact that the interaction 
changes from repulsive to attractive. 

It should be pointed out that the behavior of all the 
multiplets considered above is directly related to the 
proton-neutron effective interaction since their members 
are of almost pure configuration. In this context, it is 
quite interesting to study the renormalization of the ef- 
fective interaction, which we have taken into account 
through second-order diagrams in Vj ow -k- By way of il- 
lustration, we report in Tables IV and V the diagonal ma- 
trix elements of the effective interaction for the particle- 
hole T^hg/ 2 vf^ 2 and particle-particle 7r/ig/ 2 i / <?9/2 config- 
uration, respectively. In both tables we also show the 
corresponding matrix elements of Vi ow -k and < Vph >, 
the latter representing the contribution of the particle- 
hole core-polarization diagram, the so-called "bubble" . 
Here, we do not consider the contributions coming from 
either other second-order diagrams or folded diagrams, 
both of which we have found irrelevant for the present 
discussion. From Table IV, we see that the matrix ele- 
ments of Vph, although not very large in magnitude, are 
quite relevant in determining the pattern of the consid- 
ered particle- hole multiplet. In fact, they are repulsive for 
states with minimum and maximum J and attractive for 
J = J max — I. It is worth noting that the repulsive con- 
tribution to the Jmax state is essential to place it above 
the state with J = J m i n + 1. For the particle-particle 
case, Table V shows that the core-polarization contribu- 
tion is essentially repulsive, except for the J 71 — l~ state. 
This goes in the right direction, as it reduces the spacing 
between the 1~ and 0~ states. However, as is evident 
from Fig. 2, it is not sufficient to produce the inversion 
of these two states. 

On these grounds, we have found it interesting to mod- 
ify our effective interaction for 210 Bi by simply increasing 
the diagonal matrix elements of V^h for the ithg/2vgg/2 
configuration. It turns out that all the calculated ener- 
gies go in the right direction, a factor of 2.6 being suffi- 
cient to give the correct ground state. The results of this 
calculation (represented by open squares) are shown in 
Fig. 2. It can be seen that they practically overlap the 
experimental ones. 



B. 132 Sn region: 132 Sb and 134 Sb 

The particle- hole and particle-particle nuclei 132 Sb and 
134 Sb play in the 132 Sn region the role played by 208 Bi 
and 210 Bi in the 208 Pb region. In this connection, it 
is worth mentioning that the resemblance between the 
spectroscopy of these two regions has been emphasized 
in several recent papers |27| . 

The nucleus 132 Sb has been the subject of our study of 
Ref. Q, to which we refer for details. Some multiplets 
have been also reported in Ref. together with some 
results for 134 Sb. For completeness, we report in Fig. 3 
the 7r (77/2^4^/2 an d 7i"57/2 I/ ^-^ 1 1 /2 multiplets in 132 Sb while 
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Fig. 4 shows our results for the ^57/2^/7/2 multiplet 
in 134 Sb. Note that these multiplets are just the coun- 
terparts of those considered above for 208 Bi and 210 Bi 
(aside from the doublet ng^/ivp^^). Comments on the 
behavior of the multiplets as well as on the quality of 
agreement between theory and experiment can be found 
111 Refs. HIH. In the latter paper, the diagonal matrix 
elements of V e e, Vjow-k, and V p h for the Trg 7 / 2 vd^ 2 and 
^97 /2 v fi 72 configurations are reported. We emphasize 
here that the particle-hole and particle-particle matrix 
elements of the effective interaction in the 132 Sn region 
show the same features as those evidenced above for the 
208 Pb region. In particular, we have found that also in 
this case the core-polarization plays a significant role, but 
it is not sufficient to give the right spacing between the 0~ 
and 1" states in the particle-particle multiplet of 134 Sb 
(see Fig. 4) . We have then modified the effective interac- 
tion by multiplying the diagonal matrix elements of V p h 
for the 7r<?7/2 i/ /7/2 configuration by the same factor used 
for 210 Bi. The corresponding results are represented by 
open squares in Fig. 4, where we see that they are in 
excellent agreement with experiment. 



] Sn region: 98 Ag and 



2 In 



In the 1U0 Sn region the counterpart of 2U8 Bi and 3 Sb 
is 100 In, which until now has not been accessible to spec- 
troscopic studies. We have therefore considered the two 
neighboring odd-odd isotopes 98 Ag and 102 In , focusing 
attention on the Trg^^s^ multiplet for which some ex- 
perimental information is available. The members of the 
calculated multiplet in both nuclei have been identified 
as those dominated by this configuration with the two 
remaining valence nucleons forming a zero-coupled pair. 
In Figs. 5 and 6 we report the results of our calculations 
for Ag and 102 In, respectively, and compare them with 
the experimental data (Tjj. We see that the agreement 
between experiment and theory is of the same quality as 
that obtained in the 208 Pb and 132 Sn regions, the largest 
discrepancy being 130 keV for the 5 + state in 98 Ag. The 
pattern of the calculated multiplets turns out to be simi- 
lar to that of the multiplets in 208 Bi and 132 Sb. However, 
we find that the percentage of components other than 
those we have considered to characterize the multiplet is 
rather large. For instance, it reaches about 60% for the 
7 + state in 102 In. A detailed discussion of the structure 
of the calculated states will be included in a forthcoming 



publication. 

The diagonal matrix elements of the effective interac- 
tion for the proton hole-neutron particle con- 
figuration are reported in Table VI, where we see that 
their behavior is quite similar to that occurring in the 
heavier-mass regions. 

IV. SUMMARY AND CONCLUSIONS 

We have presented here some results of a shell-model 
study of odd-odd nuclei in the close vicinity to doubly 
magic 208 Pb, 132 Sn and 100 Sn, where use has been made 
of realistic effective interactions derived from the CD- 
Bonn NN potential. 

The main aim of this work has been to gain insight 
into the effects of the proton-neutron interaction on the 
pattern of the particle-hole and particle-particle multi- 
plets. To this end, we have made a detailed comparison 
between the calculated results and the available experi- 
mental data as well as an analysis of the proton-neutron 
matrix elements of our effective interaction. It has turned 
out that the core polarization is essential for the effective 
interaction to produce a good description of all consid- 
ered nuclei. In this context, we have found out that some 
discrepancies for the particle-particle multiplets in both 
210 Bi and 134 Sb are removed when the core polarization 
contribution is increased by a factor of 2.6. This is likely 
to be traced to the coupling of the single-particle motion 
to the octupole excitation in 208 Pb and 132 Sn. 

A relevant outcome of our calculations is that for all 
the particle- hole multiplets the J mul and J max states have 
the highest excitation energy, while the state with next 
to the highest J max is the lowest, in agreement with the 
predictions of the Brennan-Bernstein coupling rule |29| . 

In summary, we may conclude that the results of the 
present study show that realistic shell-model calculations 
are able to describe with quantitative accuracy the effects 
of the proton-neutron interaction around closed shells, 
which gives confidence in their predictive power. This 
is of particular value for nuclei in the 132 Sn and 100 Sn 
regions, for which it is of utmost importance to gain more 
experimental information. 
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FIG. 1: Proton particle-neutron hole multiplets in 208 Bi. The 
calculated results are represented by open circles while the 
experimental data by solid triangles. The lines are drawn to 
connect the points. 



FIG. 2: Proton-neutron Tvhg/ 2 i / gg/2 multiplet in 210 Bi. The 
calculated results are represented by open circles and squares 
(see text for comments) while the experimental data by solid 
triangles. The lines are drawn to connect the points. 



FIG. 3: Same as Fig. 1, but for 132 Sb. 



FIG. 4: Same as Fig. 2, but for the ng 7 /2^fr/2 multiplet in 



134 Sb. 



FIG. 5: Proton hole-neutron particle irg 9 j 2 vds/2 multiplet in 
98 Ag. The conventions of the presentation are the same as 
those used in Fig. 1. 



FIG. 6: Same as Fig. 5, but for 102 In. 
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TABLE I: Proton single-particle and neutron single-hole and 
-particle energies (in MeV) for the 208 Pb region. 



ir(n,l,j) 


e 


v{n,l,j) 1 


e 


v(n,l,j) 


e 


Ohg/2 





2pi/ 2 





109/2 





I/7/2 


0.896 


I/5/2 


0.570 


0*11/2 


0.779 


0«13/2 


1.609 


2f>3/2 


0.898 


OjlB/2 


1.423 


1/5/2 


2.826 


0«13/2 


1.633 


24/2 


1.567 


2f>3/2 


3.119 


1/7/2 


2.340 


3Sl/2 


2.032 


2^1/2 


3.633 


0^9/2 


3.414 


lff7/2 


2.491 










20(3/2 


2.538 



TABLE IV: Diagonal matrix elements of V c s , VW-k, and V ph 
(in MeV) for the -Khg^vf^^ configuration in 20S Bi. 



./ 


< V eS > 


< VW_k > 


< v ph > 


2 


0.479 


0.512 


0.102 


3 


0.140 


0.194 


-0.008 


4 


0.069 


0.129 


-0.027 


5 


0.071 


0.091 


-0.001 


6 


-0.020 


0.044 


-0.055 


7 


0.151 


0.082 


0.087 



TABLE II: Proton single-particle and neutron single-hole and 
-particle energies (in MeV) for the Sn region. 



ir(n,l,j) 


e 




e 


v(n,l,j) 


e 


Off 7/ 2 





14/2 





I/7/2 





14/2 


0.962 


0ftll/2 


0.100 


2P3/2 


0.854 


14/2 


2.439 


2si/ 2 


0.332 


0^9/2 


1.561 




2.793 


14/2 


1.655 


2pi/ 2 


1.656 


2Si/ 2 


2.800 


Off 7/ 2 


2.434 


1/5/2 


2.055 










0*13/2 


2.694 



TABLE III: Proton single-hole and neutron single-particle en- 
ergies (in MeV) for the 100 Sn region. 





e 


v{n,l,j) 


e 


0ff9/2 





14/2 





lPl/2 


0.700 


0ff7/2 


0.200 


1P3/2 


2.100 


2si/ 2 


2.200 


0/5/2 


3.100 


14/2 


2.300 






0ftll/2 


2.700 



TABLE V: Diagonal matrix elements of V e g, Vow-k, and V ph 
(in MeV) for the nh 9 /2vg 9 /2 configuration in 210 Bi. 



J 


< V eS > 


< Vlow-k > 


< Vp h > 





-0.524 


-0.514 


0.020 


1 


-0.426 


-0.333 


-0.093 


o 
2 


a oca 
-0.250 


-0.306 


A A C O 


3 


-0.179 


-0.174 


-0.002 


4 


-0.089 


-0.153 


0.065 


5 


-0.117 


-0.138 


0.025 


6 


-0.027 


-0.084 


0.061 


7 


-0.126 


-0.157 


0.034 


8 


0.015 


-0.038 


0.059 


9 


-0.275 


-0.286 


-0.003 



TABLE VI: Diagonal matrix elements of Vj.fi, Vi ow -k, and V p h 
(in MeV) for the Tvg~^ 2 vd 5 / 2 configuration in 100 In. 



J 


< VefT > 


< Vlow-k > 


< Vph > 


2 


0.994 


0.891 


0.197 


3 


0.328 


0.335 


0.020 


4 


0.171 


0.249 


-0.038 


5 


0.196 


0.175 


0.021 


6 


-0.007 


0.112 


-0.102 


7 


0.487 


0.205 


0.244 



